A modular functor which is universal for 
quantum computation 

Michael Freedman^, Michael Larsen-l^, and Zhenghan Wang^ 

February 1, 2008 



t Microsoft Research, One Microsoft Way, michaelf@microsoft.com 
I Indiana Univ., larsen@math.indiana.edu and zhewang@indiana.edu 

Abstract 

We show that the topological modular functor from Witten-Chern- 
Simons theory is universal for quantum computation in the sense a 
quantum circuit computation can be efficiently approximated by an 
intertwining action of a braid on the functor's state space. A com- 
putational model based on Chern-Simons theory at a fifth root of 
unity is defined and shown to be polynomially equivalent to the quan- 
tum circuit model. The chief technical advance: the density of the 
irreducible sectors of the Jones representation, have topological impli- 
cations which will be considered elsewhere. 
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1 Introduction 



The quantum computer was Feynman's [Fey] last great idea. He understood 
that local "quantum gates" , the basis of his model, can efficiently simulate 
the evolution of any finite dimensional quantum system and by extension 
any renormalizable system. The details of the argument are given in [LI]. 
Topological quantum field theories (TQFTs), although possessing a finite 
dimensional Hilbert space, lack a Hamiltonian — the derivative of time evolu- 
tion on which the Feynman- Lloyd argument is based. In [FKW], we provide 
a different argument for the poly-local nature of TQFTs showing that quan- 
tum computers efficiently simulate these as well. Here we give a converse 
to this simulation result. The Feynman- Lloyd argument is reversible, so we 
may summarize the situation as: 

(1) finite dimensional quantum systems, 

(2) quantum computers (meaning the quantum circuit model QCM [D] [Y]), 

(3) certain topological modular functors (TMFs). 

Each can efficiently simulate the others. We wrote TMF above instead 
of TQFT because we use only the conformal blocks and the action of the 
mapping class groups on these — not the general morphisms associated to 
3-dimensional non-product bordisms. 

We would like to thank Alexei Kitaev for conversations on our approach. 

2 A universal quantum computer 

The strictly 2-dimensional part of a TQFT is called a topological modular 
functor (TMF). The most interesting examples of TMFs are given by the 
SU(2) Witten-Chern- Simons theory at roots of unity [Wi]. These exam- 
ples are mathematically constructed in [RT] using quantum groups (See also 
[T][Wa]). A modular functor assigns to a compact surface S (with some 
additional structures detailed below) a complex vector space K(S) and to a 
diffeomorphism of the surface (preserving structures) a linear map of V{Tj). 
In the cases considered here T^(S) always has a positive definite Hermitian 
inner product <,>h and the induced linear maps preserve <,>h, i-e. are 
unitary. The usual additional structures arc fixed paramctcrizations of each 
boundary component, a labeling of each boundary component by an clement 
of a finite label set C with an involution ^ : C ^ C, and a Lagrangian sub- 
space L of ifi(E,Q) ([T][Wa]). Since our quantum computer is built from 
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quantum-S'f/(2)-invariants of braiding, and the intersection pairing of a pla- 
nar surface is 0, L = Q) and can be ignored. The parameterization of 
boundary components can be dropped. (The essential information which en- 
hances the Kauffman bracket to the Jones polynomial is remembered by the 
"blackboard framing" of the braid.) The involution " is simply the identity 
since the S'f/(2)-theory is self-dual. In fact, we can manage by only consid- 
ering the S'?7(2)-Chern-Simons theory at g = e~,r = 5 and so our label set 
will be the symbols {0, 1, 2, 3}. Note that in our notation, labels the trivial 
representation, not 1. Since we are suppressing boundary parameterizations, 
we may work in the disk with n marked points-thought of crushed bound- 
ary components. Because we only need the "uncolored theory" to make a 
universal model, each marked point is assigned the label 1, and the bound- 
ary of the disk is assigned the label 0. We consider the action of the braid 
group B{n) which consists of diffeomorphisms of the disk which leave the n 
marked points and the boundary set-wise invariant modulo those isotopic to 
the identity. The braid group has the well-known presentation: 

B{n) = {ai, ■ ■ ■ , an-i\ aiOja^^aJ^ = id if |z - j| > 1 

aiCTjai = ajOiOj if \i - j\ = 1}, 

where ai is the half right twist of the i-th marked point about the i + 1-st 
marked point. 

To describe our fault-tolerant computational model "Chern-Simons5" CSS, 
we must deal with the usual error arising from decoherence as well as a novel 
"qubit smearing error" resulting from imbedding the computational qubits 
within a modular functor super-space. To explain our approaches we initially 
ignore all errors; in particular formula (|lD is a simplification valid only in the 
error-free context. 

The state space Sk = (C^)*^^ of our quantum computer consists of k 
qubits, that is the disjoint union of k spin=i systems which can be described 
mathematically as the tensor product of k copies of the state space of 
the basic 2-level system, = span(|0 >, |1 >). For each even integer k, we 

will choose an inclusion Sk ^ V{D'^, 3k marked points) = V{D'^,3k) and 
show how to use the action of the braid group B{3k) on the modular func- 
tor to (approximately) induce the action of any poly-local unitary operator 
\J : Sk ^ Sk- That is we will give an (in principle) efficient procedure for 
constructing a braid b = 6(U) so that 

ioV = V{b)oi. (1) 
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To see that this allows us to simulate the QCM, we need to explain: (i) what 
we mean by the hypothesis "poly-local" on U, (ii) what "efficient" means, 
[iii) what the effect of the two types of errors are on line (|I|), and (iv) what 
measurement consists of within our model. 

We begin by explaining how to map Sk into V and how to perform 1 and 
2 qubit gates. 

Let D be the unit 2-dimensional disk and 

r 11 12 13 21 22 23 lOA; + 1 lOA; + 2 10A; + 3l 

\ lOOfc' lOOfc' lOOfc' lOOfc' lOOfc' lOOA;' " ' ' lOOA; ' lOOA; ' lOOfc J 

be a subset of 3k marked points on the x-axis. Without giving formulae the 
reader should picture k disjoint sub-disks D^, 1 < i < k, each containing 
one clump of 3 marked points in its interior (these will serve as qubits) and 

further ^ 2 ^ disks -Djj, 1 < ^ < J < A;, containing Di and Dj, but with 

Dij n Di = 0,1 ^ i 01 j (which will allow 2-qubit gates). Strictly speaking, 
among the larger subdisks, we only need to consider -Di,i+i, 1 < i,i + 1 < k, 
and could choose a standard (linear) arrangement for these but there is no 
cost in the exposition to considering all Dij above which will correspond in 
the model to letting any two qubits interact. Also, curiously, we will see 
that any of the numerous topologically distinct arrangements for the {Dij} 
within D may be selected without prejudice. 

We define to be the SU{2) Hilbert space of k marked points in the 
interior with labels equal 1 and / label on dD. We need to understand the 
many ways in which arises via the "gluing axiom" ([Wa]) from smaller 
pieces. The axiom provides an isomorphism: 

V{X Y) = 

tPall consistent labclings 

iV{X,l)^V{Y,l), (2) 

where the notation suppress all labels not on the 1-manifold 7 along which X 
and Y are glued. The sum is over all labelings of the components of 7 satisfy- 
ing the conditions that matched components have equal labels. According to 
S'f/(2)-Chern-Simons theory [KL], for three-punctured spheres with bound- 
ary labels a, b, c, the Hilbert space Vabc = C if 

{i) : a + b + c = even, 
(ii) : a<b + c,b<a + b,c<a + b (triangle inequalities) (3) 
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(m) : a + 6 + c < 2(r - 2); 

and Vabc = otherwise. The gluing axiom together with the above informa- 
tion allows an inductive calculation of V^, where the superscript denotes the 
label on dD. We easily calculate that 

dimV^ = 2, dimKf = 1, dimV^" = 5, dimV^^ = 8. (4) 

Line (||) motivates taking V{Di, its 3 pts, boundary label 1) =: Vi = 
as our fundamental unit of computation, the qubit. We fix the choice of an 
arbitrary "complementary vector" v in the state space of D\ U^^l Di 

V G V{D\ U^^^ Di, all boundary labels = 1 except boundary of D = 0) 
= : Complement (To keep this space nontrivial, we have taken k even.) Using 
f , the gluing axiom defines an injection: 

K : (C^)^'^ - ®UV. ^ (®t,\/.) ® Complement "'^^ C" (5) 

This composition iy determines what we will serve as our computational 
qubits within the modular functor V^^. The reader familiar with [FKW] will 
notice that we use here a dual approach. In that paper, we imbedded the 
modular functor into a larger Hilbert space that is a tensor power; here we 
imbedded a tensor power into the modular functor. 

The action of -8(3) on Di yields 1-qubit gates, whereas two qubit gates 
will be constructed using the action of B{Q) on Di j. Supposing our quantum 
computer Sk is in state s, a given v as above determines a state iv{s) = 
s ® V E V^,.. Now suppose we wish to evolve s by a 2-qubit gate g acting 
unitarily on Cf ® and hj id on Cf , I ^ i oi j . Using gluing axiom (Q) and 
the inclusion (5), we may write 

s = ^U®M/i, (6) 

h 

where {th} is a basis or partial basis for ® Cj and Uh G (8>/^j,jC^, so 
s <^ V = J2h(^h ® Uh) ® V. Decomposing along 7 = dDij, we may write 

V = ao <^ Po + a2 <S) P2, where G V [Di j\{Di U Dj), e on 7) , e = or 2 and 
/3, G V{D\{Ui^ijDi U Dij), e on 7, and on dD). Thus 

s^v = '^th^Uh®ao® l3o + '^th®Uh®a2<S) P2, (7) 

h h 
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An element of -8(6) applied to the 6 marked points in Di U Dj C Dij 
acts via a representation p° ® =: p on V^{Dij,Q pts) © V^{Dij,6 pts), 
where the superscript denotes the label appearing when the surface is cut 
along 7. In particular -8(6) acts on each factor th ® ao and th ® «2 in 
(0). Note th (S> ao belongs to the summand of V'^{Dij,6 pts) corresponding 
to boundary labels on d(^Dij\{Di U -D^)) = 0, 1, 1. There is an additional 1- 
dimensional summand corresponding to boundary labels 0,3, 3- with 0,1,3 and 

0. 3.1 excluded by the triangle inequality (ii) in (^ above. Similarly th <8) 02 
belongs to the summand of V'^{Dij, 6 pts) with boundary labels=2,l,l. There 
are additional summands corresponding to (2,1,3), and (2,3,1) of dimensions 
2 each. 

Ideally we would find a braid b = b{g) G -8(6) so that p^{h){th ® cto) = 
gth ® ao and p^{h){th ® 012) = gth ® a2. Then referring to (^ we easily check 
that 

p(6) is(g)v) = J2 iiah) ® Uh) (S) V, (8) 

h 

1. e. p{b) implements the gate g on the state space Sk of our quantum com- 
puter. In practice there are two issues: (i) we cannot control the phase of the 
output of either p° or p^, and (ii) these outputs will be only approximations 
of the desired gate g. The phase issue (i) leads to a change of the compli- 
mentary vector f — s> f ' as follows as seen on line (^ below. This is harmless 
since ultimately we only measure the qubits. 

si^v = ^ t/i ® M/i © ao /9o + ^ th®Uh®a2® P2 

h h 

J| gate 

si®v = uJo ^ gth <® Uh I® ao I® Po + UJ2 ^ gth ®Uh®a2® 132 

h h 

= ^ uogth ® M/i ® ao ® /?o + ^ ^^2gth <®Uh®a2® P2 

h h 

= "^{gth ® Uh) ® {uJoao ® (3o + uj2a2 ® 132) 

h 

=:^{gth®Uh)®v' (9) 

h 

The approximation issue is addressed by Theorem 2.1 below. 
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Theorem 2.1. There is a constant C > so that for all unitary g : ® 
C| — s> ® C^, there is a braid hi of length < I in the generators (Ti and their 
inverses a~^, 1 < i < n — 1, so that: 

\\uJop\bi) - g ® idiW + \ \uj2p\bi) - g ® td^W < e(> 0) (10) 

for some unit complex numbers (phases) Ui, i = 0,2 whenever e satisfies 

1<C-Q\ (11) 

We use II II to denote the operator norms and the subscripts on id 
indicate the dimension of the orthogonal component in which we are trying 
not to act. 

The main work in proving Theorem 2.1 is to show that the closure of the 
image of the representation p : -8(6) — > U(5) x U(8) contains SU (5) x SU (8). 
Once this is accomplished the estimate ([lO|) follows with some exponent > 2 
from [Ki] and the refinement to exponent=2 which will appear in [CN] fol- 
lowing a suggestion of the first author of the present paper. Also as explained 
in [Ki] there is a poly{\) time classical algorithm which can be used to con- 
struct the approximating braid hi as a word in {cXi} and {cr~^}. The density 
theorem is the substance of Section 4. 

The action p{h) "approximately" executes the gate g on Sk but not in 
the usual sense of approximation since the state space iv{Sk) itself is only 
approximately g invariant. To convert this "smearing of qubits" to errors of 
the type considered in the fault tolerant literature, after each g is approx- 
imately executed by p{b) we measure the labels around U^^idDi to project 
the new state p{b){s ® v) into the form s' ® v,s' G Sk, with probability 
1 — 0(e^), |s' — s| < 0{e). With probability O(e^) the label measurement 
around dDi does not yield one; in this case V^{Di;3 pts) =: Vi^i^i^i = 
has collapsed to Vi, 1,1,3 = C and it is as if a qubit has been "traced out" of 
our state space. More specifically, if the label 3 is measured on dDi, we re- 
place V^{Di, its 3 marked pts) with a freshly cooled qubit V^{D', 3 pts) with 
a completely random initial state — an ancilli — which we have been saving for 
such an occasion. The reader may picture dragging Di off to the edge of the 
disk D and dragging the ancilli D'^ in as its replacement (and then renaming 
D' by Di.) The hypothesis that such ancilli are available is discussed below. 
The error model of [AB] is precisely suited to this situation; Aharanov and 
Ben-Or show in Chapter 8 that a calculation on the level of "logical" qubits 
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can be kept precisely on track with a probability > | provided the ubiquitous 
errors at the level of "physical" qubits are of norm < 0{e) (even if they are 
systematic and not random) and the large errors (in our case tracing a qubit) 
have probability also < 0{e) for some threshold constant e > 0. For this, 
and all other fault tolerant models, entropy must be kept at bay by ensuring 
a "cold" stream of ancilli |0 >'s. In the context of our model we must now 
explain both the role of measurement and ancilla. 

Given any essential simple closed curve 7 on a surface S, the gluing 
formula reads: 

V{J:) = (BiecV{J:cut,J) (12) 
so "measuring a label" means that we posit for every 7 a Hermitian operator 



with eigenvalues distinguishing the summands of the r.h.s. of (12) above. 
For a more comprehensive computational study, we would wish to posit that 
if 7 has length = L, then can be computed in poly(L) time. For the 
present purpose we only need that -ff^,7 = dDi or dDij can be computed 
in constant time. Beyond measuring labels, we hypothesize that there is 
some way of probing the quantum state of the smallest nontrivial building 
blocks in the theory. For us these are the qubits = Vi, 1,1,1 = C^. Fix a basis 
{|0 >, |1 >} for Vi, 1,1,1 and posit for each Di, I < i < k, with label 1 on its 
boundary, an observable Hermitian operator cr* : V^j. —>■ V^j^ which acts as the 

Pauli matrix ^ ^ ^ in a fixed basis {|0 >, |1 >} for that qubit. This 

is our repertoire of measurement: is used to "unsmear physical qubits" 
after each gate and the cr^'s to read out the final state (according to von 
Neumann's statistical postulate on measurement) after the computation is 
completed. 

In fault tolerant models of computation it is essential to have available 
a stream of "freshly cooled" ancilli qubits. If these are present from the 
start of the computation, even if untouched, they will decohere from errors 
in employing the identity operator. In the physical realization of a quantum 
computer unless stored zeros were extremely stable there would have to be 
some device (inherently not unitary!) for resetting ancilli to |0 >, e.g. a 
polarizing magnetic field. As a theoretical matter unbounded computation 
requires such resetting. In a topological model such as V{T,) it is not un- 
reasonable to postulate that |0 >G Vi, 1,1,1 = ^^(-0^,3 pts) is stable if not 
involved in any gates. An alternative hypothesis is that there is some mech- 
anism outside the system analogous to the polarizing magnetic field above 
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which can "refrigerate" ancilh in the state |0 > until they are to be used. We 
refer below to either of these as the "fresh ancilli" hypothesis. To correct the 
novel qubit smearing errors, we already encountered the need for ancilli in a 



random state p = { ^ i . This state, of course, is easier to maintain. 



Let us now return to line (|l]). Let U be the theoretical output of a 
quantum circuit C of (i.e. composition of ) gates to be executed on the 
physical qubit level so as to fault-tolerantly solve a problem instance of length 
n. We assume the problem is in BQP and that the above composition has 
length < poly{n). Actually, due to error, C will output a completely positive 
trace preserving super-operator O, called a physical operator. Now simulate 
C in the modular functor V a gate at a time by a succession of braidings and 
if^-measurements. With regard to parallelism (necessary in all fault tolerant 
schemes), notice that disjoint 2 qubit gates can be performed simultaneously 
if Di ^nDi/ji = 0. For example this can always be arranged in the linear QCM 
for gates acting in -Dj^+i and -Djj+i provided i + + and i ^ j, 

and even this model is shown to be fault tolerant [AB]. As noted above, 
the complementary vector v G Kompioment evolves probabilistically as the 
simulation progresses . Different f 's will occur as a tensor factor in a growing 
number of probabilistically weighted terms. These new f-values are in the 
end unimportant; they simply label a computational state (to be observed 
with some probability) and are never read by the output measurements a*. 

Now the two main theorems: 

Theorem 2.2. Let QCM denote the exact quantum circuit model. Suppose 
M is a problem instance in BQP solved by a circuit C of length poly(L) where 
L is length(M). Let CSS denote the model based on the SU(2)-Chern-Simons 
modular functor of braids at the fifth root of unity e~ which we have de- 
scribed in this section: uncolored Sk-strand braids, and o"* measurements, 
and ^^fresh ancilli". The braid group acts on the modular functor and within 
the functor one may identify k-qubits Sk. These actions together with la- 
bel measurement 's define a probabilistic evolution of the initial (possibly 
mixed) state a G Sk. This evolution, defined gate-wise, evolves the mixed 
state a (S> f G of the modular functor to a new (probabilistic mixture of 
) state(s) (3. Performing a\- measurements on (3 samples from the mixture 
drawing out a state (3i = ai® vi and observing (according to von Neumann 
measurement) only the ai factor. With probability > | the observations cor- 
rectly solve the problem instance M. The number of marked points to be 
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braided (—3k) and the length of the braiding exceed the size of the original cir- 
cuit C by at most a multiplicative poly(log(L)) factor. Taken in triples, they 
represent the "physical qubits" of the [AB] fault tolerant model, thus CS5 
provides a model which efficiently and fault tolerantly simulates the compu- 
tations of QCM. We note that the use of label measurements introduces 
non-unitary steps in the middle of our simulation. 

Proof: The structure of the proof rehes heavily on Chapter 8 [AB] to 
reduce the QCM to a hnear quantum circuit (with state space S^) enjoying a 
very hberal error model (small systematic errors plus rare trace over qubit). 
In the final state /? = X] Pil^u each j3i admits a tensor decomposition according 
to the geometry: D = (UjDj) U (complement), but along the k boundary 
components UidDi all choices of labels 1 or 3 may appear. So if we write 
A — oii^vi we must remember that associated to I is an element [/] e {1,3}'^ 
which defines the subspaces in which ai and vi lie and that A lies in the 
corresponding [/] sector of the modular functor. All occurrences of the label 
3 correspond to a C tensor factor, C = V^{Di, 3 pts) C V{Di, 3 pts) whereas 
the label 1 corresponds to a factor. Thus in the [AB] framework each label 
3 corresponds to a "lost" or averaged qubit according to our replacement 
procedure < — > D' . Losing an occasional qubit from the computational 
space Sk is the price we pay to "unsmear" Sk within the modular functor. 
Theorem 2.1 implies that for a braid length = 0{^) a qubit will be lost with 
probability C(e^) and if no qubit is lost the gate will be performed with error 
0{e) on pure states. Factoring a mixed state as a probabilistic combination 
of pure states and passing the error estimate across the probabilities we see 
that the 0{e) error bound holds on the super-operator trace norm as well. 
Thus for e sufficiently small (estimated < 10~^ in Chapter 8 [AB]), observing 
(at random) ai amounts to sampling from an error prone implementation of 
the quantum circuit C. The error model is not entirely random in that the 
approximation procedure used to construct fe^ will have systematic biases. 
This implies that the 0{e) errors introduced in the functioning of each gate 
are not random and must be treated as "malicious". Fortunately the error 
model explained in Chapter 8 [AB] permits the small error to be arbitrary as 
long as the large error, e.g. qubit losses, occurs with a probability dominated 
by a small constant independent of the qubit and the computational history, 
as they do in our CSS model. This completes the proof of Theorem 2.2 
modulo the proof of the density Theorem 4.1. 

We may define a variant of our model CSS, "exact Chern-Simons at e^". 
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ECS5, in which we assume that all the braid groups act exactly (no error) on 
the modular functor V. Such a hypothesis is not outrageous since a physical 
implementation of a topological theory may itself confer fault tolerance, in 
that topological phenomena are inherently discrete. The only difference in 
the algorithm for modeling the QCM in ECS5 is the simplification that 
measurements are not performed in the middle of the simulation, but only 
at the very end, prior to reading out the qubits with measurements. 

Theorem 2.3. There is an efficient and strictly unitary simulation of QCM 
by ECS5. Thus given a problem instance M of length L in BQP, there is a 
classical poly (L) time algorithm for constructing a braid b as a word of length 
poly(L) in the generators a.^, 1 < i < poly(L)=3k. Applying b to a standard 
initial state, ipmiuai G V^{D, 3k), results in a state ip final G 3k), so that 

the results of on dD^ followed by a\ measurements on i) final correctly solve 
the problem instance M with probability > 6. 

Proof: In the quantum circuit model C for M (implied by the problem 
lying in BQP) count the number n of gates to be applied. Use line (11) to 
approximate each gate g hj a braid b of length / so that the operator norm 
error \\p{b) — g\ \ of the approximating gate will be less than The com- 
position of n braids which gate-wise simulate the quantum circuit introduces 
an error on operator norm < 0.1. It follows that our two measurement steps 
will return an answer (nearly) as reliable as the original quantum circuit C: 

projects to V^{D, 3 pts) with (more than) 90% probability and the sub- 
sequent probabilities of (t\ measuring |0 > or |1 > differ from C by less than 
10%. 

Remark: Theorem 2.2 and 2.3 are complementary. One provided addi- 
tional fault tolerance — fault tolerance beyond what might be inherent in a 
topological model — but at the cost of introducing intermediate non-unitary 
steps (i.e. measurements). The other eschews intermediate measurements 
by and so gives a strictly unitary simulation, but cannot confer additional 
fault tolerance. It is an interesting open problem whether fault tolerance and 
strictly unitary can be combined in a universal model of computation based 
on topological modular functors. 
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3 Jones' representation of the braid groups 



A TMF gives a family of representations of the braid groups and mapping 
class groups. In this section, we identify the representations of the braid 
groups from the SU(2) modular functor at primitive roots of unity with the 
irreducible sectors of the representation discovered by Jones whose weighted 
trace gives the Jones polynomial of the closure link of the braid [J1][J2]. To 
prove universality of the modular functor for quantum computation, we only 
use this portion of the TMF. Therefore, we will focus on these representa- 
tions. 

First let us describe the Jones representation of the braid groTips explicitly 
following [We] . To do so, we need first to describe the representation of the 
Temperley-Lieb- Jones algebras A^ji. Fix some integer r > 3 and q = . 

k -k 

Let [k] be the quantum integer defined as [k] = . Note that [—k] — 

q^Z —q-7- 

-[k], and [2] = g^ + g^. Then p := [2]^ = q+q + 2 = Acos'^{^). The algebras 
Afs^n are the finite dimensional C*— algebras generated by 1 and projectors 
Ci, • • • , e„_i such that 

1. ef = Cj, and e* = e^, 

2. eiei±iei = P~^ei, 

3. CiCj = CjCi if |i — j| > 2, 

and there exists a positive trace tr : U^iA/3,„ — > C such that tr{xen) — 
I3~'^tr{x) for all x G Ap.^. 

The Jones representation of Ap^^ is the representation corresponding to 
the G.N.S construction with respect to the above trace. An important feature 
of the Jones representation is that it splits as a direct sum of irreducible 
representations indexed by some 2-row Young diagrams, which we will refer 
to as sectors. A Young diagram A = [Ai, ■ ■ ■ Ag], Ai > A2 > • • ■ > is called 
a (2,r) diagram if s < 2 (at most two rows) and Ai — A2 < r — 2. Let 
/\n'^^ denote all (2, r) diagrams with n no des. Given A G A^^'''^ let Tf '"^ 
be all standard tableaus {t} with shape A satisfying the inductive condition 
which is the analogue of {in) in (3): when n, n — 1, • • • , 2, 1 are deleted 
from t one at a time, each tableau appeared is a tableau for some (2, r) 
Young diagram. The representation of Aj^^n is a direct sum of irreducible 
representations T^f'^'^ (2) ^) Young diagrams A. The representation 
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TT^x''^^ f*^^ ^ fixed (2, r) Young diagram A is given as follows: let v!;^'^'^^ be the 
complex vector space with basis {vt, t G T^'^^}. Given a generator Cj in the 
Temperley-Lieb- Jones algebra and a standard tableau t G V^^^'^^ Suppose i 

appears in t in row ri and column Ci, i + 1 in row r2 and column C2- Denote 

K,.+i ' 

atisfy the equation at^i = j + /^^ j. 

Then we define 



by dt^i = ci - C2 - (ri - ra), a*,, = ^p^f , and (3t^i = ^/ at^i{l - at^i). They 
are both non- negative real numbers and satisfy the equation at^i = j + Pi 



nf'''\ei){vt) = at,iVt + f3t,iVg^{t), (13) 

where gi{t) is the tableau obtained from t by switching i and i + 1 if gi{t) 
is in T^'^\ If gi{t) is not in T^'^\ then a^^j is or 1 given by its defining 
formula. This can occur in several cases. It follows that vr^^'^"* with respect 
to the basis {vt} is a matrix consisting of only of 2 x 2 and 1x1 blocks. 
Furthermore, the 1x1 blocks are either or 1, and the 2x2 blocks are 



Pt,i 1 - at,: 



(14) 



The identity at^i = a^^+P^^ implies that (|1^) is a projector. So all eigenvalues 
of Cj are either or 1. 

The Jones representation of the braid groups is defined by 

Pl3,ni^i) = g - (1 + q)ei. (15) 

Combining (p!5| ) with the above representation of the Temperley-Lieb- Jones 
algebra, we get Jones' representation of the braid groups, denoted still by 

where the dimension Npn = y^AcAt^,'-) dimK^^''^'' grows asymptotically as 
When |g| = 1, as we have seen already, Jones' representation pj^^n is 

unitary. To verify that p{(7i)p*{ai) = 1, note p*{cri) = q — {1 + q)e*. So we 

have p((Ti)p*(o-i) = qq + {1 + q){l + q)eie* - {1 + q)ei - {1 + q)e* = 1. We use 

the fact e* = Cj and = Cj to cancel out the last 3 terms. 

From the definition, p^^„ also splits as a direct sum of representations 

over (2, r)-Young diagrams. A sector corresponding to a particular Young 

diagram A will be denoted by PA,/3,n- 

Now we collect some properties about the Jones representation of the 

braid groups into the following: 
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Theorem 3.1. (i) For each (2, r)-Young diagram X, the representation px^p^n 
is irreducible. 

(a) The matrices p\^g,n{o'i) fori = 1,2 generate an infinite subgroup of 
U(2) modulo center for r 7^ 3, 4, 6, 10. 

(Hi) Each matrix p\,i3,n{o'i) , 1 < i < n — 1, has exactly two distinct eigen- 
values — l,q. 

(iv) For the (2,5)-Young diagram A = [4,2], n = 6, the two eigenvalues 
— l,g of every pA,/3,6(c"i) have multiplicity of 3 and 5, respectively. 

The proofs of (i) and (ii) arc in [J2]. For (iii), first note tliat the matrix 
P\,fi,n{(^i) is a diagonal matrix with respect to the basis {vt} with only two 
distinct eigenvalues — l,g. Now (iii) follows from the fact that all braid 
generators cxj are conjugate to each other. For (iv), simply check the explicit 
matrix for Pa,/3,6(ci) at the end of this section. 

Now we identify the sectors of the Jones representation with the represen- 
tations of the braid groups coming from the SU{2) Chern-Simons modular 
functor. The SU{2) Chern-Simons modular functor CSr of level r has been 
constructed several times in the literature (for example, [RT] [T] [Wa] [G]). Our 
construction of the modular functor CSr is based on skein theory [KL]. The 
key ingredient is the substitute of Jones- Wenzl idempotents for the intertwin- 
ers of the irreducible representations of quantum groups [RT][T][Wa]. This 
is the same SU (2) modular functor as constructed using quantum groups in 
[RT] (see [T] ) which is regarded as a mathematical realization of the Witten- 
Chern-Simons theory. All formulae we need for skein theory are summarized 
in Chapter 9 of [KL] with appropriate admissible conditions. Fix an integer 
r > 3. Let A — -\/^ • e~^, and s = A^, and q = A^. (Note the confusion 
caused by notations. The q in [KL] is which is our s here. But in Jones' 
representation of the braid groups [J2], q is A^. In all formulae in [KL], q 
should be interpreted as s in our notation.) The label set C of the modular 
functor CSr will be {0, 1, • • • , r — 2} and the involution is the identity. We 
arc interested in a unitary modular functor and the one in [G] is not unitary. 
We claim that if we follow the same construction of [G] using our choice 
of A and endow all state spaces of the modular functor with the following 
Hermitian inner product, the resulting modular functor CSr is unitary. 

Given a surface E, a pants decomposition of E determines a basis of ^(E): 
each basis clement is a tensor product of the basis elements of the constituent 
pants. The desired inner products arc determined by axiom (2.14) [Wa] if 
we specify an inner product on each space Vahc- Our choice of A makes all 
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constants S{a) appearing in the axiom (2.14) [Wa] positive. Consequently, 
positive definite Hermitian inner products on all spaces Vabc determine a pos- 
itive definite Hermitian inner product on The vector space Vabc of the 
three punctured sphere Pabc is defined to be the skein space of the disk Dabc 
enclosed by the seams of the punctured sphere Pabc- The numbering of the 
three punctures induces a numbering of the three boundary "points" of the 
disk Dabc labeled by {a, b, c}. Suppose t is a tangle on Dabc in the skein space 
of Dabc, and let i be the tangle on Dabc obtained by reflecting the disk Dabc 
through the first boundary point and the origin. Then the inner product 
<, >h' Vabc X Vabc — > C is as foUows: given two tangles s and t on Dabc, their 
product < s,t >h is the Kauffman bracket evaluation of the resulting diagram 
on obtained by gluing the two disks with s and i on them respectively, 
along their common boundaries with matching numberings. Extending <, >h 
on the skein space of Dabc linearly in the first coordinate and conjugate lin- 
early in the second coordinate, we obtain a positive definite Hermitian inner 
product on Vabc- It is also true that the mapping class groupoid actions in 
the basic data respect this Hermitian product, and the fusion and scattering 
matrices F and S also preserve this product. So CSr is indeed a unitary 
modular functor. 

This modular functor CSr defines representations of the central exten- 
sion of the mapping class groups of labeled extended surfaces, in particular 
for n-punctured disks D^^ with all interior punctures labeled 1 and bound- 
ary labeled m. If m 7^ 1, then the mapping class group is the braid group 
Bn- If m = 1, then the mapping class group is the spherical braid group 
SBn+i = M{0,n + 1). Recall that we suppress the issues of framing and 
central extension as they are inessential in our discussion. Also the represen- 
tation of the mapping class groups coming from CSr will be denoted simply 
by pr. 

Theorem 3.2. Let D^^ be as above. 

(1) : If m + n is even, and m ^ 1, then pr is equivalent to the ir- 
reducible sector of the Jones representation px^p^n for the Young diagram 
X = [m±^^ rn^] up to phase. 

(2) : If n is odd, and m — 1, then the composition of Pr with the natural 
map L : Bn ^ SBn+i is equivalent to the irreducible sector of the Jones 
representation p\,p,n for the Young diagram A = [^^, up to phase. 

The equivalence of these two representations was first established in a 
non- unitary version [Fu]. A computational proof of this theorem can be 
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obtained following [Fu] . So we will be content with giving some examples for 
r = 5. 

For the (2,5) Young diagram A = [2, 1], n = 3 with an appropriate 
ordering of the basis: 

' -1 
q 

?\/[3] 



P[2,l],/3,3(0-lJ 



P[2,l],/3,3((72) 



9+1 



9+1 
1 

■<?+l 



, where quantum [3] — q + q-\- 1. 



For the (2,5) Young diagram A = [3, 3], n = 6, the representation is 
5-dimensional. With an appropriate ordering of the basis, we have: 
/ -1 \ 

-1 

q. 

1^M 



P[3,3],/3,6(c^l) 



V 



P[3,3],/3,6(0-2j 



9 
9+1 
9V^ 

9+1 



9+1 
1 

■9+1 



Ja/[3] 



v 



9+1 
9\/[3] 
9+1 



9+1 
1 

"9+1 



^ / 

For the (2,5) Young diagram A = [4, 2], n = 6, the representation is 
8-dimensional. Here the inductive condition on basis elements make one 
standard tableau illegal, so the representation is not 9-dimensional as it would 
be if r > 5. This is the restriction analogous to {%%%) in (3) for the modular 
functor. With an appropriate ordering of the basis: 
/ -1 \ 

-1 

-1 



P[4,2l,/3,6(c^l) 
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4 A Density theorem 



In this section, we prove the density theorem. 

Theorem 4.1. Let p := pp^sj © p[4^2] : Bq U(5) x U(8) be the Jones 
representation of Bq at the 5-th root of unity q = e^. Then the closure of 
the image of p{Bq) in U(5) x U(8) contains SU{5) x SU{8). 

By Theorem 3.2, this is the same representation p := (B p'^ : Bq ^ 
U(5) X U(8) in the SU{2) Chcrn-Simons modular functor at the 5-th root 
of unity used in Section 2 to build a universal quantum computer. In the 
following, a key fact used is that the image matrix of each braid generator 
under the Jones representation has exactly two eigenvalues {—l,q} whose 
ratio is not ±1. This strong restriction allows us to identify both the closed 
image and its representation. 

Proof: First it suffices to show that the images of p[3^3] and p[4 2] contain 
SU{5) and SU{S), respectively. Supposing so, ii K = p{Bq) n {SU{5) x 
SU{8)), then the two projections pi : X — > SU{5) and p2 : K ^ SU{8) are 
both surjective. Let N2 (respectively Ni) be the kernel of pi (respectively P2). 
Then Ni (respectively A^2) can be identified as a normal subgroup of SU{5) 
(respectively SU{8)). By Goursat's Lemma (page 54, [La]), the image of K 
in SU{5)/Ni X SU{8)/N2 is the graph of some isomorphism SU{5)/Ni = 
SU{8)/N2. As the only nontrivial normal subgroups of SU{n) are finite 
groups, this is possible only if A^i = SU{5) and N2 — SU{8). Therefore, 
K = SU{5) X SU{8). 

The proofs of the density for p[3,3] and p[4^2] are similar. So we prove both 
cases at the same time and give separate argument for the more complicated 
case p[4^2] when necessary. 

Let G be the closure of the image of p[3_3] (or P[4,2]) in U(5) (or U(8)) 
which we will try to identify. By Theorem 3.1, G is a compact subgroup 
of U(m) (m = 5 or 8) of positive dimension. Denote by V the induced 
m-dimensional faithful, irreducible complex representation of G. The repre- 
sentation V is faithful since G is a subgroup of U(m). Let H be the identity 
component of G. What we actually show is that the derived group of H, 
Der{H) — [H,H], is actually SU{m). We will divide the proof into several 
steps. 

Claim 1: The restriction of V to H is an isotypic representation, i.e. a 
direct sum of several copies of a single irreducible representation of H. 
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Proof: As G is compact, V = (BpVp, where P runs through some irre- 
ducible representations of H, and Vp is the direct sum of all the copies of 
P contained in V. Since if is a normal subgroup, and the braid generators 
(Ti topologically generate G, the Ui's permute transitively the isotypic com- 
ponents Vp [CR, Section 49]. If there is more than 1 such component, then 
some (jj acts nontrivially, so it must permute these blocks. 

Now we need a linear algebra lemma: 

Lemma 4.2. Suppose W is a vector space with a direct sum decomposition 

W = (B^=iWi, and there is a linear automorphism T such that T : Wi Wj+i 
1 < i < n cyclically. Then the product of any eigenvalue of T with any n-th 
root of unity is still an eigenvalue of T . 

Proof: Choose a basis of W consisting of bases of IVj, i = 1, 2, ■ ■ ■ , n. If 
k is not a multiple of n, then trT'^ = 0, as all diagonal entries are with 
respect to the above basis. Let {Aj} be all eigenvalues of T. ( They may 
repeat.) Consider all values of trT"* = Yli\^ ("^ — 1)2,---) which are 
sums of m-th powers of all eigenvalues of T. These sums of m-th powers 
of {Aj} are invariant if we simultaneously multiply all the eigenvalues {Aj} 
by an n-th root of unity u: ^ (a;Ai)"' = ^a;™Ai™ = a;"^^Ai™ which is 
equal to trT"* = ^ A^"* because when m is not a multiple of n, they are 
both 0, and when m is, uj'^ — 1. These values trT"^ uniquely determine the 
eigenvalues of T, and therefore the set of the eigenvalues of T is invariant 
under multiplication by any n-th root of unity. 

Back to claim 1, if there is more than one isotypic component, then some 
Ui will have an orbit of length at least 2. It is impossible to have an orbit of 
length 3 or more by the above lemma as this will lead to at least 3 eigenvalues. 
If the orbit is of length 2 and as p(crj) has only two eigenvalues {a, by 
the lemma, {—a, —6} are also eigenvalues. It follows that a = —b which is 
impossible when q ^ —1. 

Claim 2: The restriction oiV to H is an irreducible representation. 

Proof: By claim 1, V\h has only one isotypic component. If V\h is 
reducible, then the isotypic component is a tensor product Vi (S> V2, where 
Vi is the irreducible representation of H in the isotypic component and V2 
is a trivial representation of H with dimV2 > 2. If Vi is 1-dimensional, 
then p{ai),i = 1, 2 generate a finite subgroup of U(m) modulo center which 
is excluded by Theorem 3.1. So we have diml^ > 2. Now we recall a 
fact in representation theory: a representation of a group p : G ^ GL{V) 
is irreducible if and only if the image p{G) of G generates the full matrix 
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algebra End(y). As Vi is an irreducible representation of H, the image 
p{H) generates End(Vi) Cg) id2, where the subscript of id indicate the tensor 
factor. As the elements (Xj normalize H, they also normalize the subalgebra 
End(Vi) ® id2 in End(yL ® V2). Consequently they act as automorphisms 
of the full matrix algebra End(Vi). Any automorphism of a full matrix 
algebra is a conjugation by a matrix, so the braid generators CTj act via 
conjugation (up to a scalar multiple) as invertible matrices in End(Vi) (E) id2 
modulo its centralizer. It is not hard to see the centralizer of End(Vi) (g) id2 
in End(Vi ® V2) is idi ® End(V2). Therefore, the braid generators CTj act via 
conjugation as invertible matrices in End(yi) (8) End(V2), i.e. they preserve 
the tensor decomposition. This is impossible by the following eigenvalue 
analysis. Consider a braid generator cij, its image p{ai) is a tensor product 
of two matrices each of sizes at least 2. Since p{(Ti) has only two eigenvalues, 
neither factor matrix can have 3 or more eigenvalues. If both factor matrices 
have two eigenvalues, the fact that p(crj) has 2 eigenvalues in all implies 
that the ratio of these two eigenvalues is ±1 which is forbidden. If one 
factor matrix is trivial, then p(o'j) acts trivially on this factor. As all braid 
generators arc conjugate to each other, so the whole group G will act trivially 
on this factor which implies that V is a reducible representation of G. This 
case cannot happen either, as y is an irreducible representation of G. 

Claim 3: The derived group, Der{H) — [H,H], of if is a semi-simple 
Lie group, and the further restriction of V to Der{H) is still irreducible. 

Proof: By claim 2, V\h is a faithful, irreducible representation, so H is 
a reductive Lie group [V, Theorem 3.16.3]. It follows that the derived group 
of H is semi-simple. It also follows that the derived group and the center of 
H generate H. By Schur's lemma, the center act by scalars. So V\Der{H) is 
still irreducible. 

Claim 4: Every outer automorphism of Der{H) has order 1, 2, or 3. 

First we recall a simple fact in representation theory. If V is an irreducible 
representation of a product group Gi x G2, then V splits as an outer tensor 
product of irreducible representations of Gi, i — 1,2. The restriction of V to 
Gi has only one isotypic component, and the restriction of V to G2 lies in 
the centralizer of the image of Gi. So the representation splits. 

Proof: It suffices to prove the same statement for the universal covering 
Der^^{H) of Der{H), as the automorphism group of Der{H) is a subgroup 
of the automorphism group of Der'^{H). 

For the 5-dimensional case: as 5 is a prime, Der^'^{H) is a simple group. 
It is well-known that any outer automorphism of a simple Lie group is of 
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order 1, 2, or 3. 

For the 8- dimensional case, if Der'"'^\H) is a simple group, it can be 
handled as above, so we need only to consider the split cases. If Der^'^{H) 
splits into two simple factors, then one factor must be SU{2): of all simply 
connected simple Lie groups, only SU{2) has a 2-dimensional irreducible 
representation. So the outer automorphism group is either Z2 when both 
factors are SU (2), or the same as the outer automorphism group of the other 
simple factor. Our claim holds. If there are three simple factors, they must 
all be SU{2). The outer automorphism group is the permutation group on 
three letters S^. Again our claim is true. 

Claim 5: For each braid generator cTj, we can choose a corresponding 
element cf, lying in the derived group Der{H) which also has exactly two 
eigenvalues, whose ratio is not ±1. The multiplicity of each eigenvalue of 
di is the same as that of cTj. (The choice of cxj is not unique, but its two 
eigenvalues have ratio q.) 

Proof: Since Der{H) is still a normal subgroup of G, and the braid 
generators cij normalize Der{H), so they determine outer-automorphisms of 
Der{H). By claim 4, an outer- automorphism of Der{H) is of order 1, 2, or 
3. Hence af acts as an inner automorphism of Der{H). By Schur's lemma, 
each af is the product of an element in Der{H) with a scalar, though the 
decomposition is not unique. Fix a choice for an element (7j in Der{H) . Then 
it has exactly two desired eigenvalues. 

To complete the proof of Theorem 4.1, we summarize our situation: we 
have a nontrivial semi-simple group Der^'^{H) with an irreducible unitary 
representation. Furthermore, it has a special element x whose image under 
the representation has exactly two distinct eigenvalues whose ratio is not ±1. 

For the 5-dimensional case, Der^'^(H) is a simple Lie group. Going 
through the list [MP] of pairs {G,zu), where G is a simply connected Lie 
group and zu a dominant weight. The only possible 5-dimensional irreducible 
representations are as follows: rank=l, {SU{2),Azui), rank=2, (5*^(4), ^^2), 
and rank=4, {SU{5),zUi),i — 1,4. By examining the possible eigenvalues, 
we can exclude the first two cases as follows: for the first case, suppose a, (3 
are the two eigenvalues of the above element x in SU{2), then under the 
representation Azui the eigenvalues of the image of x are a'^P^,i + j = 4, 
where i and j both are non-negative integers. The only possibility is two 
eigenvalues whose ratio is ±1. For the second case, since 5 is an odd number, 
any element in the image has a real eigenvalue. Other eigenvalues come in 
mutually reciprocal pairs. Again the only possibility is two eigenvalues whose 
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ratio is ±1. Therefore, the only possible pair is the third case which gives 
Der^'^{H) = SU{5). As V is a faithful representation of Der{H), the image 
of Der{H) is the same as that of Der^'^{H) which is SU{5). 

The 8-dimensional case for p[4 2] is similar. By [MP], we see the possible 
pairs for simply connected simple groups are (^SU (2), 7wij , [SU (3), W1+W2) , 
{^Spiniy),^-^) , (5*^(8), -oji) , (5*^272(8), ccj), z = 1,3,4 and (5'f/(8), Wj) , z = 
1, 7, where Wi is the fundamental weight. The same eigenvalue analysis will 
exclude all but the (S'?7(8), -ojj) case. The proof follows the same pattern 
as above with the following novelties. Case 2 is the adjoint representation 
of S'?7(3), if the special element x G S'?7(3) has eigenvalues {a,/?, 7}, the 
image matrix of x will have eigenvalue 1 with multiplicity 2 and all six pair- 
wise ratios of {a,/?, 7}, so they are ±1. For case 4, recall that if A is an 
eigenvalue of a symplectic matrix, so is with the same multiplicity, thus 
there are candidates for the special element x, but all such elements have the 
property that the multiplicity for both eigenvalues is 4. Notice by Theorem 
3.1 (iv), the multiplicity of the two distinct eigenvalue in p(crj) is 3 and 5, 
respectively. Case 5 is done just clS C3iSG 4. This excludes all the unwanted 
simple groups. We have to consider also the product cases. For product of 
two or three simple factors, the same analysis of eigenvalues as at the end of 
the proof of claim 2 excludes them. Actually, there are only four cases here: 
SU{2) X SU{2), SU{2) X SU{A), SU{2) x Spi^A) and SU{2) x SU{2) x SU{2). 
This completes the proof of our density theorem. 
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